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Abstract 

Contrary to common perceptions about systems produced in Single-Pomeron- 
Exchange (SPE) pp interactions, the hard diffractive process discovered at the 
CERN SppS-Colhder leads to dominant central production of Higgs bosons at the 
LHC. The rate for SPE production of Higgs bosons is calculated to be 7-9% of 
the total inclusive Higgs rate. In addition, an SPE measurement program of dijet 
events is outlined for the early days of LHC running which should answer many 
fundamental questions about the Pomeron structure and its effective flux factor in 
the proton. 
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1 Introduction 



A primary physics goal of the Large Hadron CoUider project (LHC) is the discovery 
and study of the Higgs boson, whose production is expected to be mainly via the gluon- 
fusion process, gg — > Higgs. The data analysis strategy for such interactions has long been 
an important topic in planning for the major experiments, ATLAS and CMS. Although 
certain types of background are irreducible, for example gg bb, it nonetheless seems 
self evident that a general reduction of hadronic activity in the final state must improve 
the capabihty for Higgs detection and isolation. 

In a typical interaction in which a 120 GeV Higgs boson is produced at the LHC, a 
gluon participating in the fusion process has only about 1% of the 7 TeV beam momentum. 
Thus, about 99% of the beam energy resides in other partons which are responsible for 
the overall hadronic activity. An alternate approach to inclusive Higgs study is to focus 
on the Pomeron-Exchange class of events (also known as "diffraction"), where much of 
the excess beam energy is carried away by a single proton on one or both sides of the 
interaction. Recently, much attention has been paid to Higgs production via Double- 
Pomeron- Exchange [1, 2] In this paper, because the cross section is much larger, we 
discuss the Single-Pomeron-Exchange process (SPE) where only a single proton emerges. 

p + p ^ p + X ^ p+ {HY), (1) 

The final state proton and the {H Y) system each recoil with more than 90% of the initial 
7 TeV beam momentum. Previous discussions of this reaction are in [3, 4, 5] with cross 
section predictions ranging from 1% to 25% of the total Higgs cross section. 

We now understand that diffractive processes occur because of the existence of colorless 
gluon-dominant clusters (for example, di-gluons) in the partonic sea of a proton which 
interact with the other beam proton. The existence of objects in a proton's sea with most 
likely momentum fraction of their host proton near zero is linked to the observations 
starting about 40 years ago [6] of inclusively-measured final-state protons from inelastic 
interactions which possessed nearly all the initial beam momentum. 

There is solid empirical evidence for the existence of such clusters. The UA8 Ex- 
periment [7] at the CERN SppS-CoUider observed jet production in this class of events, 
demonstrating the hard partonic structiirc of the clusters. The HI experiment [8] at the 
DESY HERA ep-Collidcr studied 7* interactions with rapidity gaps and demonstrated 
that the clusters have dominant gluon structure. 

In the present paper, we show that the hard Pomeron-proton collision which creates 
the (HY) system in React. 1 results in the Higgs emerging in the central region of 
the pp centcr-of-mass. To calculate the cross section for React. 1, wc assume effective 
factorization between the probabilities to find a Pomeron and its interaction with the 
other proton [9]: 

^ = Fjp/,{t,0 ■ <Jipp^hy{Mx) (2) 

where F]p/p{t, ^) is the Pomeron flux factor in the proton, which depends on the squared 
4-momentum transfer, t, to the proton and ^ = 1 — Xp, the Pomeron's momentum fraction 
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of its host proton. (7]Pp-^hy{.Mx) is the Higgs cross section, where Mx is the invariant- 
mass of the Pomeron-proton (or HY) system. To good approximation, M'jr = C,s, where s 
is the total squared cm. energy. We discuss the hmitations of the cross section calculation 
and outline how a series of measurements can be made at the LHC to test the fundamental 
aspects of Pomeron-exchange reactions. 

Section 2 contains a discussion of the kinematics of React. 1 and shows what can be 
learned from the data without prior knowledge of the Pomeron flux factor and its struc- 
ture function. Section 3 contains the calculation of the Pomeron-proton cross section in 
Eq. 2 and explains how the Pomeron structure function can be determined in hadronic 
interactions. Section 4 introduces the Pomeron flux factor and how it is determined em- 
pirically. Section 5 then uses it to calculate the total cross section for React. 1. Discussion 
and conclusions are in Sect. 6. 



2 Kinematics 

In the pp ccntcr-of-mass (the laboratory system), the mass of the Higgs, Mh, and 
its Feynman-x, x^, can be expressed in terms of the two interacting gluons' momentum 
fractions of their host protons, and ^2, respectively: 

=665 (3) 

< = 6 - 6, (4) 

where s — (14 TeV)^ at the LHC. If 6 is assumed to arise from an intermediate di-gluon 
state (the Pomeron), then: 

6 = PC (5) 

where ^ = 1 — is the Pomeron's momentum fraction of its host proton and f3 is the 
gluon's momentum fraction in the Pomeron. Eqs. 3, 4 and 5 are combined to give: 

xf ^ - H (6) 

Figure 1 displays Eq. 6 for React. 1. We see that, for fixed Mx (and therefore also ^), 
the longitudinal Higgs momentum in the laboratory is uniquely related to a value of /3. 
In other words, the observed distribution of P\\ reflects the Pomeron structure function. 
The precise relationship also depends on the gluon structure function of the proton, as 
discussed in the next section. To sum over all possible Mx values, knowledge of the 
Pomeron flux factor in the proton is needed, as discussed in Sect. 4. 

We also note that, for fixed Mx and for the maximum value, f3 = 1, P\\ has a maximum 
value in the laboratory. Fig. 1 also shows that the centrality of Higgs production can be 
influenced by an appropriate selection of Xp. Furthermore, for P|| = 0, and if we require 
^ < 0.1 (in order to have relatively clean Pomeron exchange), there is a minimum value 
of P that is sampled; it is /3 ~ 0.1 for Mh — 120 GeV at the LHC and increases either 
with larger mass or smaller s. 
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Figure 1: Higgs SPE kinematics at LHC for Mh =120 and 300 GeV,: Longitudinal Higgs 
momentum in the laboratory vs. observed proton Feynman-Xp for three values of a gluon's 
momentum fraction, (3, of its host Pomeron. Positive Higgs ' momentum corresponds to 
the direction of the observed final state proton. 

3 Higgs cross section in Pomeron— proton collisions 

The basis of all cross sections given in this paper is the cross section for gluon-gluon 
fusion into Higgs boson at a gluon-gluon mass-squared, s, which can be written as: 

agg^H{MH,s) = Kho-(Tlo-[s-S{s-M^)] (7) 

where the leading order Higgs cross section is [10]: 

Aq is the matrix element for the gg ^ H process, dominantly via a g-quark loop in leading 
order. The last factor in Eq. 7 is the zero- width approximation for the Higgs boson which, 
for Mh < 400 GeV, is known to be a good approximation. The factor, Khq, accounts 
for higher-order perturbative QCD corrections and is taken in the next-to-next-to-leading 
order [11]: Khq — K^nlo — 2 (3) for LHC (Tevatron) energy. 

In the Pomeron-proton center-of-mass system with fixed invariant mass, Mx, we write 
the cross section for inclusive Higgs boson production in terms of the standard gluon 
distribution function of the proton [12] taken at the scale = M]j, the gluon distribution 
function of the Pomeron (see below) and the cross section for gg ^ H using Eq. 7: 

{Mh,Mx) = Cg-J^ Fg/p{Q-Fg/jp{P)-agg^H{MH,s) d(5dip (9) 
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where: 

s ^ P^pM], (10) 

is the squared invariant mass of the 2-gluon system and (3 and are the gluon momentum 
fractions of their host Pomeron and proton, respectively. Cg — 0.8 is the estimated fraction 
of the Pomeron momentum that resides in the gluons [8] . 

Equation 10 and the (5-function in Eq. 7 allow Eq. 9 to be simplified to: 

aiPp^H{MH, Mx) = Cg ■ Kho ■ olo ■ j' Fg/p{^) ■ Fg/jp{(5) ■ ^ dp (11) 

where: r = {Mh/MxY- 

The Feynman-x of a produced Higgs boson in the Pomeron-proton center-of-mass is 
given in terms of the gluon momentum fractions as: 

XH ^ P - ip ^ P - t/P (12) 

Equation 12 allows us to write the differential cross section in xh as: 

Xh 



daiPp-.H{MH,Mx) T T 
^ = Cg . Kho ■ a,o ■ Fg,p{-) ■ Fg^P) ■ ^ 



1 + 



[xl + 4ry/^ 



(13) 



and where = 0.5 



XH + {xl + 4Tf' 



We evaluate these equations for several possible Pomeron structure functions [13]. 
Figure 2, reproduced from the original UA8 dijet paper [7], shows their experimental 
distribution of Feynman-x for the dijet system in the Pomeron-proton center-of-mass, 
"x(2-jet)". Almost all dijet systems travel in the Pomeron hemisphere^. As pointed 
out in [13], such an observation means that the Pomeron structure is much harder than 
the proton structure. The two curves in Fig. 2(a) show the shapes of the expected x{2- 
jet) distributions (with arbitrary normalization) for two examples of possible Pomeron 
structure functions: 

"^o/f: PFg/jp{P) ^6{l-py (14) 

''Hard" : PFg/]p{P) = 6 P (1 ~ PY (15) 

For both curves is Fig. 2(a), PYTHIA [14] was used to calculate hadronization effects and 
initial-state and final-state radiation; this was followed by a full simulation of the UA2 
calorimeter [15] which detected the jets with Et > S GeV. This method of comparing 
"smeared theoretical predictions" with observations is a reliable way to compare theory 
with experiment (see also comments in [16]). About 30% of the observed a;(2-jet) dis- 
tribution is harder than the hard function of Eq. 15. UA8 therefore called the Pomeron 
structure "Super-Hard" [7]. To parametrize this component, they assumed that the entire 
momentum of the Pomeron participated in the hard scattering as a gluon. In that case, 
p = 1 and ,r(2 jet) = 1 is the scattered parton distribution before hadronization and 

^There is no major acceptance loss in the first part of the proton hemiphere. 
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Figure 2: UA8 results [7] in the Pomeron-proton center- of-mass: (a) Observed 
Feynman-xp distribution of dijet systems. The positive x-axis is in the Pomeron hemi- 
sphere; (b) Results of X (2- jet) calculation in PYTHIA, assuming the entire momentum of 
the Pomeron participates in the hard scattering as a gluon. See discussion in teoct. 

detector effects; shown as the sohd curve in Fig. 2(b)^. The dotted curve shows the effects 
of hadronization and full detector simulation. When added to the solid (Hard) curve in 
Fig. 2(a), the dotted curve accounts nicely for the excess of events at large x(2-jet) in the 
experimental distribution. 

We therefore assume a Pomeron structure function to be the following weighted sum of 
a gluonic Hard structure, and a smeared 5-function^. For the latter, we use the function 
suggested by Alvero et al. [18]: 



'Hard + SuperHard" : pFg/ip{p) = 0.7 



6/3(1-/?) 



+ 0.3 



19.8/3^(1-/?) 



0.3 



(16) 



The essential property of this type of function required by the UA8, HI and ZEUS data 
[7, 8, 19] is that it extends to P — 1. The details of its shape do not greatly influence 
the Higgs predictions presented here. For all three types of Pomeron structure, their 
integrals are unity; in other words, a momentum sum rule is assumed to be satisfied 
and the constituents of the Pomeron carry all its momentum. The momentum sum rule 
is discussed further in Sect. 6. The evaluation of Eq. 13 using the Pomeron structure 
functions in Eqs. 14, 15 and 16 leads to the curves shown in Fig. 3 for fixed Mx — 3 TeV 
and Mh = 120 GeV. 

Following the arguments of [13], plots of experimental distributions as in Figs. 2 and 3 
directly compare the relative hardness of the Pomeron structure with the proton structure. 



^We note that because the di-gluon mass is not unique for jet formation, the sohd curve in Fig. 2(b) 
displays the features of the gluon component of the proton structure function 
^For discussion of a i5-function component in the Pomeron, see [17]. 



5 







1 1 1 1 1 
Hard + SuperHard 


Hard 






- Soft 










Dserved proton 


Mx = 3 TeV 






1 Mn = 120GeV 













-1 -0.5 0.5 1 

Higgs x^^ in IPp c. m. 



Figure 3: Differential cross section vs. Feynman-xp of a 120 GeV mass Higgs boson in 
the Pomeron-proton center-of-mass in React. 1, calculated assuming Mx = 3 TeV and 
the three Pomeron structure functions given in Eqs. 14, 15 and 16. Positive xh is in the 
observed proton direction. 



With more common final states at the LHC, such as tt, bb and dijets it will in practice be 
possible to perform high statistics determination of the Pomeron structure function as a 
function of Q^. 

Of course, following the arguments of Sect. 2, it is possible to also obtain the Pomeron 
structure function from the longitudinal momentum distributions in the laboratory using 
data at fixed Xp = 1 — ^ values. The results from Eq. 13 and Fig. 3 can be transformed 
to the laboratory frame using the Lorentz transformation: 

This leads to the distributions of longitudinal Higgs momentum in the laboratory shown 
in Fig. 4. Figure 4 also explicitly shows that, by selecting recoil proton momentum ranges 
using the Roman pots, the mean Higgs longitudinal momentum can be "tuned" . 



4 The Pomeron flux-factor 

In order to relate our calculated aiPp-,HY to observed cross sections, the probability 
to find a Pomeron in the proton at fixed ^ and t must be included. In this section we 



6 



> 0.02 

0) 

CD 

£0.01 



T3 



0.01 


0.01 



0.01 



M^^ = 120 GeV 



B, = 0.005 (Mx = 1 TeV) 



^ = 0.020 (Mx = 2TeV) 



4 = 0.046 (Mx = 3TeV) 



Hard 4 = 0.082 (Mx = 4 TeV) 

— Hard + SuperHard 



Ql 4—-'" 1 \ 

-400 -200 200 400 

Higgs P,, in pp c. m. (GeV) 



Figure 4: Higgs momentum along the beam axis in the laboratory frame (positive values 
are in the observed proton direction) for four values of Mx and ^. Mh = 120 GeV. The 
upper limit in Higgs momentum for each set of curves corresponds to (5 — 1, as seen in 
Fig. 1. 

summarize our present empirical knowledge of the Pomeron flux factor and, in the process, 
also point out its limitations and where new measurements will help. 

The Pomeron flux factor has been determined for hadronic interactions by using the 
inclusive inelastic diffractive process: 

p + p p + X (18) 

and its corresponding pp reaction. Such analyses are based on the assumption that the 

observed cross sections are the product of two factors, the Pomeron flux factor and the 
relevant Pomeron-proton total cross section, o'jpp(Mx), at invariant mass, Mx- For Re- 
act. 18, the measured cross section is written as: 

= Fw/,{t,0 ■ ajp,^x{Mx). (19) 

Refs. [21, 22] contain the results of fltting Eq. 19 to all available data on React. 18 from 
the CERN Intersecting Storage Rings [23] to the CERN SppS-Collider [21, 24], covering 
the energy range, s/s = 23 GeV to 630 GeV and over the extended ranges, |t| < 2 GeV^ 
and C < 0.09. 

The self-consistent set of fits reported in [21, 22] are equivalently important as those 
made to the diffractive DIS events by the HI [8] and ZEUS [20] collaborations. The fits 
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show that factorization between Pomeron flux and Pomeron cross section, as embodied 
in Eq. 19, describes all available data at low C,- This is all the more remarkable, since the 
empirical flux factors differ in ep and pp interactions because of their different effective 
Pomeron trajectory intercepts [9]; 

All fits in [21, 22] use the following standard Regge form for the fiux factor [25] and 
the Pomeron-proton total cross section^. Except for the Pomeron Regge trajectory which 
we discuss below, the fitted parameters in [21] are shown: 

Fjp/p{t,0 = K ■ . e(i-i±°-2)* . e^-Mt) (20) 



Ka'^lis') = (0.72 ± 0.10) ■ (4.0 ±0.6)(s')-°-''] mbGeV-^. (21) 

where s' = Mx- Equation 19 shows that, in performing fits to data, the scale factor, 
K, of the flux factor can not be separated from a scale factor, uo, of the cross section. 
Therefore, the value of the fitted product, Kgq, is given in Eq. 21. With the exponents^ 
of s' — fixed at 0.10 and -0.32, respectively, Ka^^^i^s') requires the presence of both 
Pomeron exchange and Reggeon exchange terms. With |Fi(t)p in Eq. 20 set equal to the 
Donnachie-Landshoff [25] form factor^, the additional exponential factor is required. 

Finally, concerning a(t), Ref. [22] demonstrated that, contrary to Regge theory, in 
which Regge poles are fixed, the data on React. 18 demonstrate that a{t) is an "effective 
trajectory". That is, it changes with interaction energy. Its intercept at t = decreases 
and its slope flattens as energy increases. As predicted by Kaidalov et al. [28], the reason 
for this is likely the result of complications due to multi-Pomeron-exchange effects. In 
this connection, we note that the highest energy at which a fitted effective trajectory has 
been obtained [22] is at ^/s — 630 GeV: 

a{t) = 1 + e + a't + a"t^ = 1.035 + 0.165 1 + 0.059 (22) 

It is shown in [22] that this effective trajectory is also compatible with the pubhshed 
function which was fit to the CDF data [29] at ^/s = 2 TeV. 

Since we have no information as to the effective values of a{t) at the LHC energy, all 
cross section calculations reported in this paper are given for two values of the intercept, 
q;(0) = 1.035 and 1.000. The latter value was suggested by Schuler & Sjostrand [30] for 
use at ultra-high energies. Since all evidence shows that the effective Pomeron trajectory 
intercept decreases with increasing energy, it seems like a reasonable choice to make. In 
any case, our results do not depend greatly on which of these two values is used. 

Figure 5(a) shows the .^-dependence of the flux factor of Eq. 20 integrated over all t. 
The choice of the K value assumed is perhaps the weakest aspect of the predictions, for 
there is neither a measurement nor a reliable prediction for K. The best we can do is 
use tlie \'alue of K wlvicli Doiiiiachie & Landshoff [31] extracted from elastic-scattering 

^Eq. 21 is the same function typically used to fit the s dependence of other total cross sections [26, 27]. 
^In this formula, "s"' stands for "s'/sq" , where sq = 1 GeV^. 

-^iW - 4m|-t ■ (l-t/0.7l)2 
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Figure 5: Dependence on ^ and Mx of: (a) Flux factor in Eq. 20 integrated over all t, 
using the Pomeron effective trajectory in Eq. 22 with e — 0.035, and then 0.0. K = 0.78 
GeV^^ is assumed, (b) The Pomeron-proton to Higgs cross section, calculated for the 
two indicated Higgs masses using Eq. 11 and the Hard + SuperHard Pomeron structure 
function, (c) Cross section for React. 1 (the product of the solid curve in (a) and each 
curve in (h)), multiplied by a factor of 2 to account for final state proton detection in both 
arms. 

data, namely K = 0.78 GeV^^, although given the observed s-dependence of the effective 
Pomeron trajectory discussed above, we see no reason why this value of K should also 
apply precisely to inelastic diffraction where multi-Pomeron-exchange effects may be 
different. We therefore regard it as an approximation. However, we discuss tests of this 
issue in Sect. 6. 



5 Total Higgs cross section via SPE 

The total Higgs cross section via the SPE mechanism can now be calculated. Fig. 5(b) 
shows aiPp^niMn, Mx) plotted vs. Mx and ^ for two values of Mh- The product of the 
flux factor and aiPp^niMn, Mx) is the t-integrated version of Eq. 2 and is shown in 
Fig. 5(c), multiphed by a factor of 2 to account for the use of Roman pot spectrometers 
in both arms. aip{0) — 1.0 is assumed for both parts (b) and (c) of Fig. 5. Finally, Fig. 6 
shows the ^-integral of the differential cross section in Fig. 5 vs. Higgs mass. At LHC, 
SPE Higgs production (with recoil protons measured in either arm) is between 7% and 
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Figure 6: Upper curve is the total inclusive Higgs boson cross section vs. Higgs mass, 
calculated for the gluon-gluon fusion process. Central curves are the LHC cross section of 
Fig. 5(c) integrated over ^ < 0.1; for the two indicated values ofa]p{0). The lower curves 
are the same, but for the Tevatron. 

9% of the total Higgs cross section, depending on Mh- 

The Higgs longitudinal momentum in the laboratory integrated over all Mx is shown in 
Fig. 7. This figure gives us global views of the centrality of Higgs production integrated 
over all ^ in the range: < ^ < 0.1. The distributions do not depend greatly on the 
detailed shape of the Pomeron structure function at large (3. 




6 Discussion Conclusions 

We have drawn attention to the method used by the UA8 Collaboration [7, 13] to 
determine the Pomeron structure function. The longitudinal momentum distribution of 
the system produced via hard scattering refiects the relative hardnes of the Pomeron and 
proton structure functions. The analysis can be caried out either in the Pomeron proton 
cm. or in the laboratory although, in the later case, it must be done at fixed value of 
Xp = 1 — ^. Progress in further understanding the Pomeron will be possible at the LHC by 
making high-statistics measurements of more common final states than those with Higgs 
bosons, for example dijets, tt or bb. 

One of the parameters most poorly understood is the normalization, K, of the Pomeron 
fiux factor in the proton. To further understand this issue, the following ratio can be 
measured over a large range of t and ^ : 

^ _ ^'^pp-*p{dijet+Y) ^23^ 
A(7pp_>p(x) 
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Figure 7: Differential cross sections at LHC and Tevatron for React. 1 vs. P\\ of the Higgs 
in the laboratory frame for the two indicated Higgs masses and two assumed Pomeron 
structure functions. 

If factorization between flux factor and IPp cross sections remains valid, the ratio, 7^, will 
be found to be independent of t at fixed ^. If that is the case, then the flux factor cancels 
out of the ratio and we have: 

Numerator and denominator of Eq. 24 have been multiplied by K, because the denom- 
inator is then a measurable quantity (by fitting Eqs. 20, 21 and 22 to the data). The 
numerator of the right-hand side of Eq. 24 can then be extracted and compared with the 
calculated dijet cross section. 

This comparison allows K to be determined although, strictly speaking, a possible 
violation of the momentum sum rule in the calculated dijet cross section means that the 
quantity measured is actually "fK'\ where / < 1.0 if the sum rule is violated. UA8 
carried out this procedure [32] with their limited statistics dijet data. For a pure gluonic 
Pomeron, they found: 

fK = 0.38 ±0.13 (25) 

We have already pointed out that the Donnachie-Landshoff value, K — 0.78 GeV~^ 
which was extracted from clastic scattering data may not be precisely applicable to in- 
elastic diffractive data, because of different multi-Pomeron-exchange effects. Thus, if the 
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momentum sum rule is valid, Eq. 25 is a measure of K for inelastic diffraction. On the 
other hand, it could reflect a product of the two effects. In either case, taking Eq. 25 at 
face value would mean that all cross sections presented in this paper arc overestimated 
by a factor of 2.0 ± 0.7. Studying fK obtained in this way for different final states and 
over a wide range of kinematics should shed additional light on its correct interpretation. 

In conclusion, we have shown that Hard Diffraction Scattering with Single-Pomeron- 
Exchange leads to Higgs production dominantly in the central region. The extent to which 
the decrease in hadronic activity characteristic of such reactions will actually improve 
Higgs identification and isolation is best shown with the aid of full event simulation using 
the PYTHIA software package [33]. Such studies are underway. We have also briefly 
outlined an initial program for SPE production of dijet events which should lead to major 
improvements in our knowledge of the Pomeron structure and its flux factor in the proton. 
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